Astronomy & Astrophysics manuscript no. Bonomo CoRoTlOb 


©ESO 2010 


June 16, 2010 





o 

(N 

1—5 

in 



Oh 
I 

O ■ 



Transiting exoplanets from tlie CoRoT space mission* 

XI. CoRoT-10b: a giant planet in a 13.24 day eccentric orbit 
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ABSTRACT 

Context. The space telescope CoRoT searches for transiting extrasolar planets by continuously monitoring the optical flux of thousands of stars in 
several fields of view. 

Aims. We report the discovery of CoRoT-lOb, a giant planet on a highly eccentric orbit (e = 0.53 ± 0.04) revolving in 13.24 days around a faint 
(V= 15.22) metal-rich KIV star 

Methods. We use CoRoT photometry, radial velocity observations taken with the HARPS spectrograph, and UVES spectra of the parent star to 
derive the orbital, stellar and planetary parameters. 

Results. We derive a radius of the planet of 0.97 ± 0.07 Rj„p and a mass of 2.75 ± 0.16 Mj„p. The bulk density, pp = 3.70 + 0.83 g cm"\ is ~ 2.8 
that of Jupiter The core of CoRoT-lOb could contain up to 240 of heavy elements. 

Moving along its eccentric orbit, the planet experiences a 10.6-fold variation in insolation. Owing to the long circularisation time, Tci,c > 7 Gyr, a 
resonant perturber is not required to excite and maintain the high eccentricity of CoRoT-lOb. 

Key words, planetary systems - stars: fundamental parameters - techniques: photometric - techniques: spectroscopic - techniques: radial veloc- 
ities. 
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1. Introduction 

CoRoT is the pioneer space mission dedicate d to the detec tion 
of extrasolar planets via the transit method (iBaehnl 120031: see 
lAuvergne et al.l l2009l for a detailed description of the instru- 
ment and its performance). To date it has led to the discov- 
ery of ten extrasolar planets: CoRoT-7b, the first super-Earth 
with m easured radius and mass dLeger et a l.l'2009f: ^Ouel oz et alJ 
l2008l) : three inflated hot J upiters CoRoT-lb (Barae et al. l|2008h 
CoRoT-2b (" Alonso et all 12008a') and CoRoT-5b (Rau eretalJ 
L2009); CoRoT-3b, "th e first secure inhabitant of the brown dwarf 
desert" dPeleuil et alJ l2008): two Jupiter-like planet s with an or- 
bital period of approxima tely 9 days, CoRoT-4b ( Aigr ain et alj 
l2008HMoutou et al.ll2008h and CoRoT-6b (Fridlund et al.i 2010h:_ 
the l ong-period temperate giant planet CoRoT-9b (DeegetaD 
2010), and the hot Super-Neptune CoRoT-8b ( Horde et al] 
20ld). 



as HD 147056 b (ahas HAT-P-2b e ^ 0.52 and orbital period 
P = 5.63 davs: iBakos et alJl2007t iPal et aUlIoTol). H D 17156b 



Here we report the discovery of the giant planet CoRoT- 
10b that orbits its parent star in 13.24 days, moving along a 
highly eccentric orbit with e = 0.53 + 0.04. It is therefore 
one of the few known transiting planets with e > 0.5, such 



Send offprint requests to: A. S. Bonomo 
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* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribution 
of Austria, Belgium, Brazil, ESA (RSSD and Science Programme), 
Germany and Spain. 



(e = 0.67 and P = 21.22 days; 'Barbieri et al."2007', '2009) and 
HD 80606b (e ^ 0.93 and P ^ 1 1 1 .44 days ; Moutou et al. 200§ 
iWinn et al.l2009l:fHebrard et al.l2010l) . Eccentric transiting plan- 
ets include also those with a smaller eccentricity, notably the 
recendy discovere d WASP-8b (e - 0.31 and P - 8.16 days; 
|Queloz etal.''2010) and HAT-P-15b (e = 0.19 and P = 10.9 
days; Kovacs et al. 2010 ), the massiv e planet X 0-3b (e - 0.26 
and P = 3.19 davs; IJohns-Krull et al . 2008; Winn et al.l 120081) 
and the two Neptunes GJ 436b (e = 0.15 and P - 2.64 d ays; 
iGillon et al.ll20 07a b; A lonso et al.ll20"0 8b; Bean et al l2008[) and 
HAT-P-1 lb (e = 0.20 and P = 4.89 davs; lBakos et a0201ot) . 



Transiting planets in eccentric orbits are very intriguing 
and interesting objects as they allow us to study ongoing tidal 
dissipation and its impact on the planet radius (Ibgui et al] 
2010), atmospheri c circulation in the case o f a strong varia- 
tion in insolation (iLangton & LaughTinll2007h . and the dynam- 
ical orbital evolution including the gravitational interaction be- 
tween planets in a multiple system (planet-planet scattering; 
e.g., Marzari & Weidenschilling 2002) or the secular influence 
of a p ossible distant stellar companion (Kozai mechanism; iKozail 
Il962h . 
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Fig. 1. The quiescent light curve of CoRoT-10 binned at 512 s 
showing ten transits of the giant planet CoRoT-lOb. Jumps due 
to hot pixels were removed by means of an iterative 3-sigma 
clipping. 

Table 1. CoRoT-10 IDs, coordinates and magnitudes. 



CoRoT window ID 


LRc01.E2_1802 




CoRoT ID 


100725706 




USNO-A2 ID 


0900-14919216 




2MASS ID 


19241528+0044461 




GSC2.3 ID 


NIMR021985 




Coordinates 






RA (J2000) 


19:24:15.29 




Dec (12000) 


00:44:46.11 




Magnitudes 






Filter 


Mag 


Error 


B" 


16.68 


0.14 


V 


15.22 


0.05 


r'" 


14.73 


0.03 




13.74 


0.03 


I* 


12.53 


0.02 


H" 


11.93 


0.03 


K" 


11.78 


0.02 



" Provided by Exo-Dat JPeleuil et aU2"009t) : 
* from 2MASS catalogue. 



oversampling modeQ. Fig. [T] shows the CoRoT-10 light curve 
with the nominal sampling of 512 s, filtered from a) outliers 
that are produced by proton impacts during the crossing of the 
South Atlantic Anomaly of the Earth's magnetic field by the 
satellite; and b) several jumps, with a typical duration shorter 
than 1 day, due to hot pixels . Corrections for the CCD zero 
offset, sky background. Earth's scattered light and jitter varia- 
tions were carried out by the latest version 2.1 of the CoRoT re- 
duction pipeline. Unl ike CoRoT-2 (Lanza et al. 2009), CoRoT-6 
(iFridlund et al.ll2010l) and CoRoT-7 (Lanza et aL201 0). the Hght 
curve of CoRoT-10 is relatively quiescent and does not show flux 
variations due to the presence of starspots and photospheric fac- 
ulae greater than a few mmag. It shows a long term decrease of 
~ 2.5% attributable to CCD ageing (Fig.lTJ. The rms of the nom- 
inal and oversampled photometric points is 0.0013 and 0.0046 in 
relative flux units, respectively. 

A total of 10 transits with a depth of ~ 1 .3% are visible in the 
light curve (Fig.lTJ. A zoom of one of the five 32 s oversampled 
transits is shown in Fig. |2] The transit ephemeris, reported in 
TablelH was derived from a linear fit to the measured transit mid 
points determined by a trapezoidal fitting of each transit. It gives 
an orbital period of 13.2406 + 0.0002 days. 



3. Ground-based follow-up observations 

3.1. Radial velocity observations 

We performed radial velocity (RV) obser vations of t he star 
CoRoT-10 with t he HARPS spectrograph dPepe et all F2002bl 
iMavor et al.ll2003l) at the 3.6-m ESO telescope (La Silla, Chile). 
HARPS was used with the observing mode obj_AB, without ac- 
quisition of a simultaneous Thorium lamp spectrum in order to 
monitor the Moon background light on the second fibre. The in- 
trinsic stability of this spectrograph does not require the use of 
lamp calibration spectra, the instrumental drift during one night 
being in our case always smaller than the stellar RV photon noise 
uncertainties. HARPS data were reduced with the on-line stan- 
dard pipeline and radial velocities were obtained by a weighted 
cross-correlation with a numerical spectral mask for a K5V star 
(Baranne 1996; I'ene el al. 2Q,02a.) . 



' data available at http://idoc-corot.ias.u-psud.ft71 



2. CoRoT observations 

The parent star of CoRoT-lOb, i.e. the CoRoT target 
LRc01_E2_1802, is a V = 15.22 star and has been observed 
in the stellar field pointing towards the constellation of Aquila 
during the first CoRoT long run LRcOl (Cabrera et al. 2009). 
Its magnitudes in several photometric bands and its coordinates 
are reported in Table[T] CoRoT observations of this target lasted 
for 142.07 days, from the 16* of May up to the 15* of October 
2007, and provided us with monochromatic (white channel) data 
(Auvergne et al. 2009) . 

Tra nsits by CoRoT-l Ob were first discovered in "alarm 
mode" (ISuraceet al.ll20 08'). i.e. while CoRoT observations were 
still ongoing, which permitted us to change the temporal sam- 
pling from 5 1 2 s to 32 s after HJD 2454305 . 1 1 . In total, 2 1 0248 
photometric measurements were obtained, 198752 in the 32 s 
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Fig. 2. One of the five 32 s oversampled transits of CoRoT-lOb. 
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Fig. 3. 7b;? panel: radial velocity measurements obtained and the 
Keplerian best-fit solution (solid line). Bottom panel: residuals 
from the best-fit. The open circles indicate the measurements 
affected by the moonlight after our correction. 



The first two measurements of CoRoT- 1 were made on June 
20080. Seventeen additional measurements were obtained from 
June to September 2OOS0. Seven of our nineteen measurements 
were strongly contaminated by the moonlight: the radial veloc- 
ity of the Moon was close to that of CoRoT-10 and affected both 
the RV measurements and the bisector lines. We developed a 
software correction using the Moon spectrum simultaneously ac- 
quired on fibre B: it consists in subtracting the cross-correlation 
function (CCF) of fibre B, containing the Sun spectrum (re- 
flected by the Moon), from the CCF of the fibre A, containing the 
stellar spectrum. The correction was applied when the two CCF 
peaks were close in radial velocity. For CoRoT-10, corrections in 
the range between 50 and 550 m s"' were applied for seven mea- 
surements. To be conservative, we added quadratically 30 ms ' 
of systematic errors to these seven corrected measurements. 



trie planetary companion. We fitted the data with a Keplerian 
orbit using the CoRoT ephemeris P - 13.2406 days and 
r„ = 2454273.3436 HJD (see Table Hi. The derived eccen- 
tricity and argument of periastron are e - 0.53 + 0.04 and 
a; = 218.9 + 6.4 deg. The other orbital parameters are reported 
in Table H] The standard deviation of the residuals to the fit 
cr{0 -C) = 29 m s"' is comparable to the mean RV uncertainty. 

To examine the possibility that the RV variation is due to a 
blended binary scenario - a single star with an unresolved and 
di luted eclipsi ng binary -, we followed the procedure described 
in lBouchv et a l. (2008) based on checking both the spectral line 
asymmetries and the dependences of the RV variations on dif- 
ferent cross-correlation masks. These two checks excluded that 
the RV variation is caused by a blended binary and allowed us to 
establish the planetary nature of CoRoT-lOb. The bisector varia- 
tions are shown in Fig.|5] 
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Fig. 5. Bisector variations (span of the bisector slope) as a func- 
tion of orbital phase {left pannel) and radial velocity {right pan- 
nel). Bisector error bars are estimated as twice the radial veloc- 
ity uncertainties. No bisector effect is visible for the moonlight- 
corrected measurements (open circles) indicating the good qual- 
ity of our correction. 
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Fig. 4. Phase-folded radial velocity curve of CoRoT-10 and the 
Keplerian best-fit solution (solid line). The open circles indicate 
the measurements affected by the moonlight after our correction. 

The radial velocities are listed in Table |2] and displayed in 
Fig. |3] and |4] The phase-folded radial velocity curve shows a 
variation in phase with the CoRoT transit period. It is compat- 
ible with the reflex motion of the parent star due to an eccen- 



3.2. Photometric observations 

Photometric observations during and outside the transit were 
carried out at the 1 .20 m telescope of the Observatoire de Haute 
Provence during the nights of the 18* of June and the 15* of 
July 2009, respectively. Such observations are complementary to 
the radial velocity measurements and are required to definitively 
exclude the possibility that the transits detected in the CoRoT 
light curve could be produced by a background eclipsing bi- 
nary wh ich contaminates the CoRoT aperture mask of the star 
(Deeg et al. 2009). The latter covers ~ 16" x 14" on the sky 
and contains two main background contaminants: one is 3.78 
mag fainter in V and is located at ~ 6.5" from CoRoT-10 to- 
ward South- West; the other one is 4.22 mag fainter in V and is 
~ 7.5" toward North-East (see Fig. |6ll. Differential photometry 
shows an "on-target" transit with the same depth as observed in 
the CoRoT light curve. None of the two contaminants exhibits 
significant flux variations that can mimic the transits observed in 
the CoRoT-10 Hght curve. 



4. Transit fitting 



2 //Ai?P5 program 08 l.C-0388 

3 i/Ai?P5 program 083.C-0186 



In order to perform the transit fitting, first of all we filtered the 
raw light curve from outliers due to impacts of cosmic rays. 



4 



A. S. Bonomo et al. 2010: CoRoT-lOb 



Table 2. Radial velocity measurements of CoRoT-10 obtained 
by HARPS. HJD is the Heliocentric Julian Date. The seven ex- 
posures which were affected by the moonlight are labelled with 
an asterisk and 30 m s ' was quadratically added to their errors. 



HJD 


RV 


±lo- 


BIS 


exp. time 


S/N/pix. 


-2400000 


[kms-'] 


[kms-'] 


[kms-'] 


[s] 


(at 550nm) 


54640.85815* 


15.386 


0.038 


0.021 


3600 


7.4 


54646.74608 


15.049 


0.035 


0.087 


1800 


5.1 


54986.78569* 


15.322 


0.049 


0.081 


3600 


4.5 


54989.77932* 


14.947 


0.035 


-0.048 


3600 


9.2 


54990.81302* 


14.923 


0.038 


0.004 


3600 


7.1 


54994.81853* 


15.545 


0.036 


-0.014 


3600 


8.4 


54995.80369* 


15.479 


0.032 


-0.050 


3600 


10.8 


55022.79584* 


15.503 


0.035 


0.015 


3600 


7.0 


55063.69761 


15.507 


0.032 


0.055 


3300 


5.4 


55068.67093 


15.079 


0.018 


-0.024 


3600 


7.7 


55069.68356 


14.947 


0.019 


0.006 


3300 


7.6 


55070.69173 


15.146 


0.021 


0.070 


3000 


7.4 


55071.51030 


15.414 


0.019 


-0.074 


3000 


7.8 


55072.54251 


15.456 


0.020 


-0.014 


3000 


7.5 


55073.57439 


15.492 


0.015 


-0.030 


3000 


9.1 


55074.61800 


15.508 


0.022 


-0.073 


3600 


7.7 


55075.53585 


15.468 


0.029 


-0.018 


3000 


6.4 


55077.62048 


15.398 


0.032 


0.114 


3000 


5.7 


55078.61399 


15.343 


0.024 


0.018 


3600 


7.6 



1 arcmin 




E 



Fig. 6. The sky area around CoRoT-10 (the brightest star near the 
centre). Left: R-filter image with a resolution of ~ 2.5" taken 
with the OHP 1.20 m telescope. Right: image taken by CoRoT, 
at the same scale and orientation. The jagged outline in its centre 
is the photometric aperture mask; indicated are also CoRoT's x 
and y image coor dinates and position of nearby stars from the 
Exo-Dat database (iDeleuil et al.ll2009l) . 



Based on pre-launch o bservations stored in the Exo-Dat cata- 
logue dPeleuil et a l. 2009), we estimated the flux contamination 
from the two faint background stars which fall inside the CoRoT- 
10 photometric mask to be 5.5+0.3% (see Fig.|6]l. We subtracted 
such a value from the median flux of the light curve (75559 
e-/32 s), which makes the transits slightly deeper by ~ 7 • 10 "* in 
relative flux. We then fitted a parabola to the 5 h intervals of the 
light curve before the ingress and after the egress of each transit 
in order to correct for any local variations. We disregarded two 
of the ten CoRoT-lOb transits, precisely the second transit and 
the ninth, as their shape was deformed by hot pixels. Finally, we 
folded the light curve using the ephemeris reported in Table |4] 
and binning the data points in bins of 8 • 10"'' in phase, cor- 



responding to ~ 1.5 min (Fig. |7]i. The error on each bin was 
computed as the standard error of the data points inside the bin. 

Trans i t best - fit wa s performed following the formalism of 
iGimenej (l2006i l2009l) and fixing the eccentricity and the argu- 
ment of the periastron to the values derived from the Keplerian fit 
of the radial velocity measurements (see Sect. 13. lb . The four free 
parameters of the transit model are: the transit centre; the phase 
of the end of transit egr ess 62 in the reference sys t em fo r ec- 
centric orbits defined by iGimenez & Garcia-Pelavd (Il983h and 
shown in their Fig. 1 (see also Gimenezll2009l) ; the ratio of the 
planet to stellar radii k - Rp/R,, and the inclination / between 
the orbital plane and the plane of the sky. The two non-linear 
limb-darkening coeflicients M+ = Ug + ut, and m_ = Ma-Mi,0were 
fixed for two reasons: first, the relatively low signal-to-noise ra- 
tio of the transit light curve, which does not allow us to con- 
strain either m+ or m_ within reasonable error bars; secondly, the 
degeneracy between the inclination and the two limb darkening 
coeflicients in the case of a high impact parameter b, as in our 
case b - 0.85 + 0.03 (see TableHJi. The adopted limb-darkening 
coeflicients Ua and ui, for the CoRoT bandpass were determined 
following the procedure in Sing (2010). However, while the lat- 
ter takes into account only the stellar emergent intensity for the 
values of /i > 0.05, we considered all the 17 values available in 
the ATLAS model grids0down to /i = 0.01. Our choice is moti- 
vated by the fact that the transit of CoRoT-lOb is grazing, which 
implies that the variation of the specific intensity close to the 
limb of the stellar disc must be considered properly when mod- 
elling the transit shape. The derived limb-darkening quadratic 
coeflicients are Ua - 0.51 and ut, - 0.21, which give m+ = 0.72 
and M_ - 0.3. 

The best-fit tran sit paramete rs were found by using the al- 
gorithm AMOEBA (iPress et al.i ri992) and changing the initial 
values of the parameters with a Monte-Carlo method to find the 
global minimum of the x^- Our best-fit of the phase-folded and 
binned transit light curve is shown in Fig. [T] Fitted and derived 
transit parameters are listed in Table |4] together with their 1- 
si gma errors est i mated u sing the bootstrap procedure described 
in 'Alonso et alJ (I2008ah which takes also the correlated noise 
into account (cf. lAlonso et al . 2008a, Sect. 3). Uncertainties on 
the eccentricity, the argument of the periastron and the contami- 
nation were also considered for the estimation of the errors of the 
derived transit parameters b, a/Rt, a/Rp, MJ" /Rt and p», where 
b is the impact parameter, a the semi-major axis of the planetary 
orbit, Rt and M* the stellar radius and mass, Rp the planet ra- 
dius, and ft the stellar density (see Table |4]l. The fitted value of 
the transit centre is consistent with zero within 1-sigma, which 
confirms the good quality of the transit ephemeris (Sect.|2]i. 



5. Stellar and planetary parameters 

The spectral analysis of the parent star was performed with a 
high-resolution UVES spectrum acquired on the 30''' of July 
2OOS0 We used the Did mode (390+580) and a sHt width of 
0.8", achieving a resolving power of ~55 000. The total exposure 
time was 4 h leading to a signal-to-noise ratio per pixel S/N~120 
at 5500 A. 



* u„ and ui, are the coefficients of the limb-darkening quadratic law: 
= 1 -Ma(l -ii)-Ub(l —j-if', where /(I) is the specific intensity 
at the centre of the disk and yu = cos 7, 7 being the angle between the 
surface normal and the line of sight 

^ http://kurucz.harvard.edu 

' t/V£5 program ID: 083.C-0690(A) 
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Fig. 7. Top panel: Phase-folded light curve of 8 transits of CoRoT-lOb. The bin size corresponds to 1.5 min and the 1-sigma error 
bars on each bin are estimated as the standard error of the data points inside the bin. The solid line shows our best-fit transit model. 
Bottom panel: the residuals from the best-fit model. 



Table 3. Abundances of some chemical elements for the fitted 
lines in the UVES spectrum. The abundances refer to the solar 
value and the last column reports the number of lines used. 



Saturated interstellar Na D lines in the HARPS spectra indi- 
cate a significant absorption along the line of sight. Converting 



Element 


Abundance 


No. lines 


Cai 


0.21 ±0.12 


5 


Til 


0.38 ±0.11 


15 


Vi 


0.63 ±0.12 


10 


Cri 


0.34 ±0.13 


7 


Fe I 


0.26 ±0.10 


77 


Fe n 


0.26 ±0.11 


5 


Co I 


0.41 ±0.10 


7 


Nil 


0.35 ±0.10 


23 


Si I 


0.37 ±0.10 


8 



To derive the stellar atmospheric parameters, we first de- 
termined the V sin it - 2 + 0.5 km s '. To that purpose, we 
selected a few HARPS spectra that were not contaminated by 
the Moon reflected light. This series of spectra was set at rest 
and co-added. To carry out the det ailed spectral analysis, we 
made use of the VWA (lBrunttll2009i) software package and ob- 
tained: Teff = 5075 + 75 K, log 4.65 + 0.10 cm s^^ and 
[Fe/H]= 0.26 + 0.07 dex. The surface gravity value was checked 
with usual indicators: Na I D lines around 5890 A, Mg I b lines, 
and Ca I lines at 6122, 6162 and 6439 A. The abundances of 
several chemical elements are listed in Table[3] The elements for 
which we could only measure a few lines are not reported. 

The absence of noticeable emission in the core of the Ca II 
H & K lines supports the low magnetic activity of CoRoT-10 
indicated by its quiescent light curve (see Fig.[T]). 



initudes (Table 


mi in 


Bessell & Brett) 


1988") 



ing the ( J - K) colour with that expected by Kurucz mod- 
els for the CoRoT-10 spectral type and metallicity, we found 
a colour excess E{ J - K) ^ 0.24. This corresponds to an ex- 
tinction of Ay ^ 1.39 ma^ in agreement with reddening maps 
(ISchlegel etalJ [T998^. Using the Pogson formulc0, the stellar 
distance was estimated to 345 + 70 pc. 

From the transit best-fit we derived a stellar den sity of 
Pt = 3.32 + 0.70 g cm-^ i.e. 2.3 5 ± 0.50 Po. C ESAM (More 
2008) and STAREVOL dPalaciosJIprivate comm unication: Sies 
2006) models of stellar evolution do not foresee any evolu- 
tionary track that matches the above-mentioned stellar density, 
given the effective temperature and the metallicity of CoRoT- 
10. Specifically, they predict an upper-limit of 2.55 g cm"^ 
(1.79 po), compatible at 1. 1-sigma with the stellar density de- 
rived from the transit fitting. The mass and radius of the star 
provided by the CESAM evolutionary tracks are respectively 
equal to M. = 0.89 + 0.05 Mg and = 0.79 ± 0.05 Rq. The 
corresponding surface gravity, logg=4.59 + 0.06 cm s"^, is in 
good agreement with the spectroscopic value. The stellar age 
constraints are relatively weak but favor values smaller than 3 
Gyr. 



' AvIEiJ -K) = 5.82 ± 0.1 ICo3l2000h 

V - My = 5 logd - 5 + Ay, where My was determined from the 
bolometric magnitude, given the T^s and R, of CoRoT-10, and the BC 
taken from |http://kurucz.harvard.edu for the atmospheric parameters of 
CoRoT-10 
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From the aforementioned stellar radius and mass, we deter- 
mined the radius of the planet Rj, = 0.97 + 0.07 Rjup and its 
mass Mo = 2.75 + 0.16 Mjup. The bulk density, pp = 3.70 + 



0.83 g cm ^, is ~ 2.8 that of Jupiter. 



6. Discussion and conclusions 

We report the discovery of CoRoT-lOb, a transiting planet on 
a highly eccentric orbit (e = 0.53 + 0.04) with a mass of 
2.75 + 0.16 Mjup and a radius of 0.97 + 0.07 Rjup. It orbits a 
metal rich ([Fe/H] = 0.26 + 0.07) KIV star with a visual mag- 
nitude V = 15.22 in 13.24 days. Fig. [8] shows the position of 
CoRoT-lOb in the eccentricity -period diagram of the known ex- 
trasolar planets and highlights its peculiarity as it belongs to the 
class of the few transiting exoplanets with highly eccentric orbits 
(e > 0.5) among which HAT-P-2b, HD 17156b and HD 80606b. 

To investigate the internal structure of CoRoT-lOb, 
we computed planetary evolution models with CERAM 
(iGuillot & Morel 1995) under the standard hypothesis that the 
planet is made up of a central rocky core of varia ble mass and o f 
an overlying envelope of solar composition (e.g.. lGuillotll2008l) . 
The results in terms of planetary size as a function of system 
age are shown in Fig.|9] The coloured regions (red, blue, yellow- 
green) indicate the constraints derived from the stellar evolution 
models at 1 , 2 and 3 cr level, respectively. Assuming a zero Bond 
albedo, we derived the equilibrium temperature of the planet 
Teq = 600 K. For this temperature, models of planet internal 
structure with a core mass of 0, 20, 60, 120, 180, 240 and 320 
Me were computed (Fig.|9l). Note that a 25% change in the equi- 
librium temperature yields a difference in the resulting planetary 
radius of less than 1%. Therefore, a Bond albedo considerably 
different from zero (up to Ab = 0.7) does not change signifi- 
cantly our results. 

CoRoT-lOb is a high density plan et with a mass an d density 
similar to those of HD 17156b (Bar bieri et alj|2009l) but with 
a higher content of heavy elements. For an age of the star and 
hence of the planetary system between 1 and 3 Gyr, CoRoT- 
10b should contain between 120 and 240 M^ of rocks in its in- 
terior (i.e. between 14 and 28% of the total mass), at 1-sigma 
level. Mixing heavy elements in the envelope rather than assum- 
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Fig. 8. Eccentricity - period diagram for the known extraso- 
lar planets (black filled symbols are transiting planets). The 
size of the symbol indicates the mass range. CoRoT-lOb, with 
P = 13.24 days and e = 0.53, is indicated by the red filled circle. 
Data from http://exoplanet.eu, 
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Fig. 9. Evolution of the radius of CoRoT-lOb (in Jupiter units) as 
a function of age, compared to constraints inferred from CoRoT 
photometry, spectroscopy, radial velocimetry and CESAM stel- 
lar evolution models. Red, blue and green areas correspond to the 
planetary radii and ages that result from stellar evolution mod- 
els matching the inferred p, - T^ff uncertainty ellipse within Icr, 
2cr and 3cr, respectively. Planetary evolution models for a planet 
with a solar-composition envelope over a central dense core of 
pure rocks of variable mass are shown as solid lines. These mod- 
els assume a total mass of 2.75 Mjup and an equilibrium temper- 
ature of 600 K (corresponding to a zero Bond albedo). They de- 
pend weakly on the assumed opacities, and uncertainties due to 
the atmospheric temperature and planetary mass are negligible. 



ing that they are all containe d in the core may y ield a reduction 
of these numbers by ~ 30% (iBaraff'e et al.ll2008 ). This number 
is uncertain however because it does not account for the increase 
in opacity in the outer radiative zone that would have the oppo- 
site effect (Guillot 2005) . In any case, CoRoT-lOb is found to be 
extremely enriched in heavy elements, sug gesting that its form a- 
tion probably required giant collisions (see Iko ma et al.ll2006h . It 
also strengthens the observed correlation between star metallic- 
ity and heavy elements in the planet ( Guillotil2()08h . 

Moving along its eccentric orbit from the periastron to the 
apoastron, CoRoT-lOb experiences a 10.6-fold variation in inso- 
lation. Moreover, since the tidal interaction with the parent star 
is strongest around periastron, the planet i s exp ected to be in 
pseudo-synchronous rotation. Eq. 42 from iHutI (Il98lh gives a 
rotation period of the planet Pp.iot = 4.25 + 0.53 days. Time- 
dependent radiative models of planetary atmospheres could be 
used to study heating variations due to the changing star-planet 
distance at different pressure levels of the planet atmosphere 
and to predict temperature inversion s caused by s t rong h eat- 
ing around periastron, as done by liro & DemingI (1201 Oh for 
HD 17156b and H D 80606b. 

Using Eq. 6 in Matsumura et al. ('2008"), we can estimate the 
circularisation time Tchc of the planetary orbit, neglecting the 
stellar damping. The circularisation time strongly depends on 
the adopted tidal quality factor for the planet Q'^ which is not 

well known. For Ju piter, 5 ■ 10"* < gp < 2 ■ 10^ dYoder & Peald 
IT98llLainev et alj]2009l) and for mos t of the giant extrasola r 
planets 10^ < QJ, < lO' (see Fig. 2 from lMatsumura et al.l2008h . 

Considering Q'^ - 10^, we find Tcirc ~ 7.4 Gyr. Owing to the 
long circularisation time, the eccentricity of CoRoT-lOb need 
not be excited and maintained by the resonant interaction with 
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another planet. Nevertheless, the eccentric orbits of extrasolar 
planets can be explained by gravitational planet-planet scattering 
(e.g., Ichatterjee et al. 2008 and references therein). If another 
massive planet survived the violent encounters between planets 
and is currently orbiting around the parent star, it could be de- 
tected by a long-term radial velocity follow-up of the parent star, 
e.g., showing a long-term drift induced by the distant compan- 
ion. If the star has a distant companion of stellar nature, the high 
eccentricity of CoRoT-lOb could be produced by Kozai oscilla- 
tions rather than planet-planet scattering (e.g.. ITakeda & Rasiol 
Distinguishing between the two scenarios would make us 
understand the dynamical evolution of the eccentric giant planet 
CoRoT-lOb. 
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Table 4. Planet and star parameters. 



Ephemeris 


Planet orbital period P [days] 


13.2406 + 0.0002 


Planetary transit epoch T,, [HJD-2400000] 


54273.3436 ± 0.0012 


Planetary transit duration rftr [h] 


2.98 + 0.06 


Planetary occultation epoch T„^^ " [HJD-2400000] 


54276.49 ± 0.41 


Planetary occultation duration docc [h] 


1 no 1 n 1 o 
z.Oo ± U. lo 


Epoch of periastron To [HJD-2400000] 


54990.85 ± 0.08 


Derived parameters from radial velocity observations 




Orbital eccentricity e 


0.53 + 0.04 


Argument of periastron oj [deg] 


218.9 ± 6.4 


Radial velocity semi-amplitude A" [ms"'] 


301 ± 10 


Systemic velocity [kms ^] 


Ij.jjU ± u.uu/ 


O-C residuals [ms"'] 


29 


Fitted and fixed transit parameters 






0.00483 + 0.00009 


Radius ratio k = R^jR, 


0.1269 ± 0.0038 


Inclination i [deg] 


88.55 ± 0.2 




-1-0 79 


K_ (fixed) 


+030 


Derived transit parameters 




alR, ' 


31.33 ±2.15 


a/Rp 


247 ± 21 


(MJM^y'^RJR^)-' 


1.33 ±0.09 


Stellar density [gem ^] 


j.jZ ± U. /U 


Impact parameter b 


0.85 ± 0.03 


Spectroscopic parameters of the star 




Effective temperature reff[K] 


5075 ± 75 


Surface gravity log g [cgs] 


4.65 ± 0.10 


Metallicity [Fe/H] [dex] 


+0.26 ± 0.07 


Stellar rotational velocity Vsin!» [kms"'] 


2.0 ± 0.5 


Spectral type 


KIV 


Stellar and planetary physical parameters 




Star mass [Mq] 


0.89 ± 0.05 


Star radius [Rq\ ' 


0.79 ± 0.05 


Planet mass Mp [Mjjip] 


2.75 ±0.16 


Planet radius R^ [Rjup] 


0.97 ± 0.07 


Planet density pp [g cm^^] 


3.70 ± 0.83 


Planet surface gravity log gp [cgs] 


3.93 ± 0.08 


Planet rotation period Pp^ot [days] ^ 


4.25 ± 0.53 


Distance of the star d [pc] 


345 ± 70 


Orbital semi-major axis a [AU] 


0.1055 ±0.0021 


Orbital distance at periastron a^^y [AU] 


0.0496 ± 0.0039 


Orbital distance at apoastron a.^po [AU] 


0.1614 ± 0.0047 


Equilibrium temperature at the averaged distance T^q [K] ^ 


600 ± 23 


Equilibrium temperature at periastron T^' [K] * 


935 ± 54 


Equilibrium temperature at apoastron T^^° [K] * 


518 ±20 



" r„„ = 7'i,+ I -(f +(1 +csc-/)-e cosoi); 

^ phase of the end of transit egress in the reference system defined bv lGimenez & Garcia-Pelavol U983I) . 

<: a/R. = '^'""^ ■ -I = ==, where V2 is the true anomaly measured from the periastron passage at the end of transit 

^_l2__^-cos2(v2+^-f).si„2/ 

egress (see lGimene 



^ from CESAM stellar evolution models. See Sect[5]for details; 

^ assuming the planet to be in a pseudo-synchronous rotation; 

^ black body equilibrium temperature for an isotropic planetary emission. 



